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In this paper we aim to study the chemical evolution of sul- 
phur in the galactic disk, using a new optimal abundance indi- 
cator: the [S i] line at 10821 A. Similar to the optimal oxygen 
indicators, the [O i] lines, the [S i] line has the virtues of being 
less sensitive to the assumed temperatures of the stars inves- 
tigated and of likely being less prone to non-LTE effects than 
other tracers. High-resolution, near-infrared spectra of the [S i] 
line are recorded using the Phoenix spectrometer on the Gemini 
South telescope. The analysis is based on ID, LTE model at- 
mospheres using a homogeneous set of stellar parameters. The 
A 10821 [S i] line is suitable for an abundance analysis of disk 
stars, and the sulphur abundances derived from it are consis- 
tent with abundances derived from other tracers. We corrobo- 
rate that, for disk stars, the trend of sulphur-to-iron ratios with 
\Q metallicity is similar to that found for other alpha elements, 
C> supporting the idea of a common nucleosynthetic origin. 
O ■ 

, Key words, stars: abundances - stars: atmospheres - stars: 
^ ■ late-type - Galaxy: disk - infrared: stars 
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Abstract. 

1. Introduction 

To extract information about the nucleosynthesis of the a ele- 
ments (e.g., Mg, Si, S, and Ca), accurate data on their abun- 
dances for stars in the galactic halo and disk are essential. 
Sulphur is one of the least studied a elements, owing to the 
paucity of appropriate spectral lines. The diagnostics tools used 
for sulphur so far are multiplets of transitions from highly 
excited levels, such as the -U6743 - 57, AA8694 - 95, and 
/L19213 - 38 A lines (corresponding to the oxygen triplet at 
7775 A). The derived abundances from these lines are sensi- 
tive to the adopted T e ff of the star and could be susceptible to 
non-LTE effects that are difficult to quantify owing to a lack of 
adequate atomic data. Thus, these multiplets may lead to un- 
certain sulphur abundances. A recent debate in the literature 
has questioned the origin of sulphur in the early universe, as 
well as whether sulphur follows the evolution of the other a 
elements or not, partly owing to these uncertainties, but also 
due to the weakness of lines used (see for example Israelian 
& Rebolo 2001; Takada-Hidai et al. 2002; Nissen et al. 2004; 
Ryde & Lambert 2004, 2005; Korn & Ryde 2005). 

Sulphur and oxygen belong to the same group in the pe- 
riodic table and have the same number of valence shell elec- 
trons, which means that their spectra show some similarities, 
and similar diagnostic lines are available for the determination 
of their elemental abundances in stars. In the same way as for 
the /16300 and /16363 A [O i ] lines from the ground state, which 
have been extensively used and discussed for the determina- 
tion of oxygen abundances in stars (see for example Asplund 
et al. 2004; Nissen et al. 2002; Garcia Perez et al. 2006), the 
corresponding [S i] lines should not be susceptible to non-LTE 
effects, nor should they be very dependent on the adopted T e g . 
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These virtues make it worthwhile to investigate the A 1 082 1 [Si] 
line as a spectroscopic tracer of stellar sulphur abundances. The 
reason why this line has not previously been used is that it lies 
beyond the reach for normal optical CCDs. The near-infrared 
region also has several other advantages for an abundance study 
of cool stars (Ryde et al. 2005). 

In this Letter we present the first attempt to determine stel- 
lar sulphur abundances based on the forbidden sulphur line at 
10821.18 A (in air). This line is an inter-combination line, an 
Ml transition between the triplet ground-state CP2) and the 
first excited singlet state ( D2). We explore its use in stars of 
the galactic disk since it is important to evaluate this new pre- 
ferred diagnostic and to scrutinize the evolution of sulphur in 
the disk based on other diagnostics. The abunda nces of sulphu r 
in disk s t ars have earlier been stud ied bv Iciegg et all fl9 8 1 ), 
iFrancoisI ( fl987l) . IChen et aD d2002l) . and lCaffau et alJ d2005t) . 

2. Observations 

The [S 1] line was observed at a high signal-to-noise ratio 
(SNR ss 200), in a sample of G and K subgiants and giants 
of the galactic disk 1 . The homogeneous set of stars was se- 
lected from the library of F5-K7 stars observed with the Elodie 
spectrometer (Soubiran et al. 1998) and are listed in Tabled 
Our high-resolution, near-infrared spectra were recorded with 
the Phoenix spectrometer mounted on the Gemini South tele- 
scope. Phoenix is a cryogenic, single order echelle spectrom- 
eter with an Aladdin 512 x 1024 element InSb detector array 
cooled to 35 K (Hinkle et al. 1998). We used a 3 pixel wide slit 
resulting in a spectral resolution of R ~ 60, 000. The spectra 
were centered on the A 10821 [S 1] line and covered the 9220- 
9260 cm -1 region. The data were obtained in service observing 
mode (Programme ID: GS-2004A-Q-13) over a period from 
5/4 2004 to 23/6 2005. 

1 For the purpose of the present paper we define all stars with 
metallicity [Fe/H] > -1 as "disk" stars. The notation [A/B] = 
log(A , A/A , B)* - log(A , A/A , B)o> where Na and are the number abun- 
dances of elements A and B, respectively. 
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Fig. 1. The spectrum of the star HD2 12943 is shown by dots. 
The full line shows a synthetic spectrum with [S/Fe] = 0.21, 
based on a model atmosphere with the following fundamental 
parameters: T eS = 4586 K, logg = 2.81 (cgs), [Fe/H]= -0.34, 
£micro = 1.5kms _1 , and^macro = 4.5kms _I . 

The data were processed in a standard manner with 
the reduction package IRAF (e.g., Tody 1993) to retrieve 
one-dimensional, continuum-normalized, and wavelength- 
calibrated stellar spectra. The raw images were trimmed and 
a normalized fiat frame was made by dividing a master fiat 
by a Legendre polynomial to fit the overall spectral shape in 
the dispersion direction. The observations were performed in 
pairs such that the star was moved approximately 5" along the 
slit between two equally long exposures. These two exposures 
were subsequently subtracted from each other, a process that 
reduces the sky background. A hot calibration star, recording 
telluric atmospheric absorption lines, was used partly for the 
wavelength calibration and partly to eliminate the telluric spec- 
tral contamination. The multiple exposures of the same star 
made it possible to apply a cosmic -ray-hit rejection algorithm 
when adding them together in the form of 1 -dimensional spec- 
tra. Due to residual fringing, a high-order polynomial was used 
to rectify and normalize the spectra. An example of a spectrum, 
with lines identified, is shown in Fig.d 

In the analysis the equivalent width of the Al 0821 [S i] line 
was measured using the deblending algorithm in IRAF, since 
the wing of a Cri line blends into the wing of the [Si] line. 
Gaussian fits were used, and we estimate the measuring un- 
certainty to be 10%, leading to an uncertainty in the derived 
sulphur abundance of the order of 0.03 dex, as estimated from 
individual measurements. The equivalent widths are given in 
Tabled 

3. Analysis 

For the analysis of the stellar spectra, ID Marcs model atmo- 
spheres (Gustafsson et al., in preparation) are computed for 
each star. Sulphur abundances were derived from the equiv- 
alent widths, using the Eqwi code, and the synthetic spectra 
were calculated with the Bsyn code. Both codes use many sub- 
routines from the Marcs code and solve the radiative transfer 



at many frequencies in each spectral line using the continu- 
ous opacity from the model atmosphere and the line opacity 
calculated from line lists. The fundamental stellar parameters 
of the observed stars, which are neede d in the calculation of 
the stellar atmospheres, are taken from ISoubiran et all £l998) 
and are also given in Table d Moreover, a microturbulence 
of 1.5kms _1 is used and the usual linearly increased (for de- 
creasing metallicity) a element enhancement for disk stars is 
assumed. The Marcs models used are hydrostatic model photo- 
spheres in spherical geometry (although not really necessary 
for the stars investigated here) and are computed assuming 
LTE, chemical equilibrium, homogeneity, and the conservation 
of the total flux (radiative plus convective, the convective flux 
being computed using the mixing length recipe). The radiation 
field used in the model generation is calculated with absorp- 
tion from atoms and molecules by opacity sampling at approx- 
imately 95 000 wavelength points over the wavelength range 
1300 A-20//m The atomic opacity sampling files are calcu- 
lated for metallicities of [Fe/H]= 0.0, -0.1, -0.25, -0.37, -0.5, 
and -0.75, and the one closest in metallicity is used. 

Molecular lines are included, but only CN (line list from 
J0rgensen & Larsson 1990, and Plez 1998, private commu- 
nications) is important for the most metal-rich stars with low 
effective temperatures. TiO is, for example, unimportant even 
for the coolest of the stars. Atomic line data are taken from the 
VALD database (Piskunov et al. 1995) with the exception of 
the transition probability of the [S i] line, which is taken from 
the NIST database version 3 (see, e.g., Martin et al. 2000). The 
wavelength and strength of the Fe i line and the position of the 
two Cr i lines in the vicinity of the [S i] line are matched to fit 
the solar lines (see Table|2}. To take additional line broadening 
into account, the synthetic spectra are convolved with a radial- 
tangential function (Gray 1992) with a macroturbulent velocity 
given in Table d as £ maC ro- The synthetic spectra are used to 
inspect the location of the continuum, the line strengths and 
widths, and line blendings. For the stars we are investigating, 
no weak atomic blends are found at the [S i] line. The sulphur 
abundances are given in Tabled 

In Fig. |2] we have plotted our derived values of [S/Fe] vs. 
[Fe/H], as well a s the the rel a tive su lphur abundances derived 
for disk stars by IChen et all d2002h . who used optical, high- 
excitation S i lines (x exc = 7.9 eV). 



4. Discussion 

We ha ve chosen to use the stars in the library of Sou biran et alJ 
dl998l) to keep to a set of stars with homogeneously determined 
stell ar parameter s. The accuracies of the parameters given in 
iKatz et all Jl998h for the stars in this library are 8T eS = 86 K, 
6logg = 0.28 dex, and 5[Fe/H] = 0.16dex. In Table|3]the un- 
certainties in the derived sulphur abundances due to uncertain- 
ties in the fundamental stellar parameters are shown. We also 
show the corresponding changes for the nearby Cr i and Fe i 
lines. Including the statistical measurement uncertainties, as 
estimated from the continuum placement and the effect of the 
SNR (both max 10% in W$), a conservative total uncertainty in 
the derived sulphur abundances is estimated to be 0.13 dex. The 
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Table 1. The stars investigated. Columns 3-6 give the atmospheric parameters of the stars, after Soub iran et alJ Jl998l) . The two 
subsequent columns provide the dates and the exposure times of the observations of the stars. The last four columns give our 
measured equivalent widths of the A 10821 [S i] line, the derived logarithmic number abundance of sulphur relative to hydrogen, 
the logarithmic sulphur abundance relative to the Sun, and finally, the logarithmic sulphur -to-me tal licity ratio relative to the solar 
values. In the second part of the table the stellar parameters for two stars from Fuhrmannl d2004l) . and the sulphur abu ndances we 
derive based on these, are presented. We have used a solar sulphur abundance of loge Q (S) = log Ns/Nn + 12 = 7.20 dChen et all 
l2002h . 



Star 


Spectral 




log(g) 


[Fe/H] f 


micro / smaero 


Date of 


Exposure time 


w s 


loge(S) 


[S/H] 


[S/Fe] 




type 


[K] 


[cgs] 




[km s" 1 ] 


observations 


M 


[mA] 








HD 3546 


G5III 


4942 


2.73 


-0.66 


1.5/5.5 


July 6, 2004 


120 


11.8 


6.79 


-0.41 


0.25 


i in i no on 




a n^T 
4LD / 


1 A1 

1.4 J 


— U.ZJ 


1.3/ J. J 


July d, ZUU4 


ou 


jb.y 


£ no 
o.yo 


— U.zZ 


U.Uj 


HD 404fSf) 


K1III 


4741 


2.00 


-0.50 


1.5/5.0 


Dec 1Q 7004 


480 


18.3 


6.70 


-0.50 


0.00 


HD 81192 


G7III 


4705 


2.50 


-0.62 


1.5/5.5 


April 6, 2004 


512 


12.7 


6.73 


-0.47 


0.15 


HD 117876 


G8III 


4782 


2.25 


-0.50 


1.5/5.0 


April 5, 2004 


480 


19.9 


6.80 


-0.40 


0.10 


HD 139195 


KOp 


4856 


2.82 


-0.04 


1.5/5.5 


May 17, 2004 


88 


13.0 


7.02 


-0.18 


-0.14 


HD 161074 


K4III 


3980 


1.73 


-0.27 


1.5/6.0 


June 22/23, 2005 


960 


60.9 


7.13 


-0.07 


0.20 


HD 168723 


K0III-IV 


4859 


3.13 


-0.19 


1.5/4.7 


May 17, 2004 


20 


12.4 


7.15 


-0.05 


0.14 


HD 184406 


K3in 


4520 


2.41 


+0.01 


1.5/5.0 


Aug. 7, 2004 


120 


25.1 


7.15 


-0.05 


-0.06 


HD 188512 


G8IVvar 


5041 


3.04 


-0.04 


1.5/5.0 


July 5, 2004 


180 


9.1 


6.98 


-0.22 


-0.18 


HD 212943 


K0III 


4586 


2.81 


-0.34 


1.5/4.5 


June 29, 2004 


84 


19.9 


7.07 


-0.13 


0.21 


HD 214567 


G8II 


5038 


2.50 


+0.03 


1.5/5.5 


June 30, 2004 


288 


19.7 


7.15 


-0.05 


-0.08 


HD 219615 


G7III 


4830 


2.57 


-0.42 


1.5/6.0 


June 29, 2004 


64 


19.6 


7.02 


-0.18 


0.24 


HD 220954 


K1III 


4664 


2.37 


-0.10 


1.5/5.5 


June 29, 2004 


48 


30.0 


7.19 


-0.01 


0.09 


HD 168723 


K0III-IV 


4921 


2.99 


-0.21 


1.09/4.7 


May 17, 2004 


20 


12.4 


7.07 


-0.13 


0.08 


HD 188512 


G8IVvar 


5110 


3.60 


-0.17 


0.92/5.0 


July 5, 2004 


180 


9.1 


7.23 


+0.03 


+0.20 



originates in the ground state, the sulphur line opacity (or the 
number of sulphur atoms in the ground state) is not very tem- 
perature sensitive. The continuous opacity of the wavelength 
region (H^) is slightly T-sensitive, but we see from Table [3] 
that the sulphur abundance derived is not very T-sensitive. Nor 
is the sulphur line opacity very sensitive to the electron pres- 
sure, although the continuous opacity is. This explains the sen- 
sitivity of the abundance derived from the observed equivalent 
width to the metallicity and especially to the surface gravity, as 
seen in Table [3] The Cr and Fe lines behave differently since 
they are from excited states, and (like H ) represent a minority 
species in the photosphere. Therefore, the derived Cr and Fe 
abundances are more temperature sensitive, but are not as sen- 
sitive to the electron pressure. It can also be noted that the [S i] 
line is not saturated, and therefore not sensitive to £ m i cr o- 

In the discussion on errors we have ignored systematic er- 
rors due to the treatment of convection in the atmospheres. A 
more proper way of treating convection is through ab initio 
3D h ydrodynamic models as described in, e.g., lAsplund et alJ 
(2000J). For the [Oi] lines. iNissen et al.l i2002h show that for 
disk metallicities the 3D effects are less than 0. 1 dex (decreas- 
ing with increasing metallicity) for their coolest dwarf (5800 
K). The main effect is the change in the continuous opacity 
and not a change in the line opacity due to cooler upper layers. 
Sulphur is expected to behave in a similar fashion. The effect is 
unknown for giants and subgiants, and we did not include it in 
the present analysis. 

In his careful analysis of nearby disk and halo stars, 
iFuhrmannl d2004l) has two stars in common with our star list, 
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Fig. 2. [S/Fe] vs. [Fe/H] plot for the stars investigated in th is pa- 
per us ing the A 10821 [S i] line (circles) and for the lChen et alJ 
(2002J) stars (crosses). Typical error bars are indicated. Squares 
indicate new sulphur abundances derived from a different set of 
stellar parameters for rj Ser and j3 Aql (see text). 



uncertainty in the metallicity is taken from lKatz et alJ 1^998), 
i.e., <5[Fe/H] =0.16 dex. 

The strengths of weak lines depend on the line and contin- 
uum opacities, xiim and^- cont , such that the equivalent width of 
a line W °c ^line/^fcont- Since sulphur is mainly in the atomic 
state throughout the photospheres and since the weak [S i] line 
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Table 2. Atomic line parameters. After the excitation energy, 
Xexc, of the lower level of the transition, which is given in eV, 
the transition probability of the line as the logarithm of the 
product of the statistical weight and the oscillator strength is 
given. 



Element 


Wavenumber 


Wavelength 


ATexc 


log Sf 






(in vacuum) 








[cm- 1 ] 


[A] 


[eV] 




Cri 


9242.264 


10819.860 


3.013 


-1.957 


Fei 


9241.082 


10821.244 


3.960 


-2.178 


[Si] 


9238.609 


10824.140 


0.000 


-8.617 


Cri 


9238.200 


10824.619 


3.013 


-1.678 



namely 77 Ser (HD168723) and Aql (HD188512). The stellar 
parameters derived by him, and the alternative sulphur abun- 
dances derived by us using his stellar parameters, are also given 
in Table [2 As we saw in our uncertainty analysis, the most 
sensitive fundamental parameter of the stellar models influ- 
encing the abundance derived from the [S 1] line is the sur- 
face gravity. In particular, the alternative surface gravity of 
HD188512 (logg = 3.04 -> 3.60) is the major cause of the 
large change in the sulphur abundance of that star (log e(S) = 
6.98 -> 7.23). Since the metallicity of HD188512 is also dif- 
ferent ([Fe/H]-0.04 — » -0.17), the position of the star in Fig. 
|5]changes quite dramatically. This again gives an indication of 
the uncertainties due to the stellar parameters only. 

From Fig. Ew e see that the agre ement with the abundance 
determination of Chen et al .1 (120021) is good given the level 
of uncertainties. IChen et al.l J2002I) showed that their S abun- 
dances correlated strongly with their Si abundances, the latter 
being a well studied a element, suggesting that the nucleosyn- 
thesis of sulphur is similar to that of the other a elements. We 
can thus corroborate that sulphur shows an evolution typical of 
an a element, and we have shown that the A 10821 [S 1] line is 
indeed a good diagnostic tool for the sulphur abundance in disk 
stars. The virtues of the line are that it is not very sensitive to 
T e g and that it should not be subject to non-LTE effects. It is 
apparent that an uncertainty in the surface gravity affects the 
abundance derived from this line. It should, however, be noted 
that the [S/Fe] ratio will be much less dependant on gravity 
if the iron abundance is determined from Fe 11 lines that react 
similarly to a change in electron pressure as the [S 1] line. The 
question whether the A 10821 [S 1] line can also be used for stars 
of the galactic halo will be addressed in a forthcoming paper. 
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